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Abstract 

Time-dependent soliton solutions are explicitly derived in a five-dimensional theory en- 
dowed with one (warped) extra-dimension. Some of the obtained geometries, everywhere well 
defined and technically regular, smoothly interpolate between two five-dimensional anti-de 
Sitter space-times for fixed value of the conformal time coordinate. Time dependent solutions 
containing both topological and non-topological sectors are also obtained. Supplementary 
degrees of freedom can be also included and, in this case, the resulting multi-soliton solutions 
may describe time-dependent kink-antikink systems. 
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1 Physical motivations 



Time-dependent soliton solutions arise typically in (1 + 1) space-time dimensions by appro- 
priately boosting static solutions. This observation is due to the pioneering work of various 
authors [1, 2, 3] (see also [4, 5] and references therein for an extended introduction) whose 
main focus was on the Sine-Gordon system. The static solutions of the (1 + 1) dimensional 
system only involved one scalar degree of freedom. Furthermore, the geometry was always 
static and Minkowskian both in the case of static and time- dependent solutions. The boosted 
solutions could be simply interpreted, from the physical point of view, as traveling solitons 
which could even illustrate, with appropriate boundary conditions (see, for instance, [3]), 
the elastic collision of two profiles. 

In (1 + 1) dimensions, static multisoliton solutions are also known to exist. The obtainable 
solutions can be categorized into two basic groups. To the first set belong solutions containing 
two (or more) topological defects. The second group of solutions contemplates two (or more) 
defects of which at least one is a non-topological soliton. In the figurative language of some 
of the authors discussing these profiles they have been dubbed as bag-like solutions [6, 7, 8]. 
According to the same language it is also possible to achieve trapping-bag configurations 
where a topological kink is trapped by a non-topological bag (see also [9, 10] where trapping 
bag solutions appear in the framework of the so-called MIT-bag and SLAC-bag models of 
extended hadrons). The interest of multisoliton solutions in (1 + 1) dimensions persists 
even in recent years and various interesting integration methods have been tailored to deal, 
specifically, with non-linear multi-field equations [11, 12, 13]. 

The main theme of the present paper will not be (1 + 1) dimensional systems but rather 
five-dimensional gravitating soliton solutions. It is nonetheless inspiring to bear in mind 
the aforementioned observations since the topic of the present paper will be time- dependent 
generalizations of static soliton solutions in five-dimensional warped space-times. 

Gravitating solitons in higher dimensions can be used to model, for instance, five- 
dimensional domain walls where the extra (bulk) coordinate is infinite and where the space- 
time geometry is consistently included as a solution of the appropriate higher-dimensional 
gravitational field equations. Single gravitating kinks in five-dimensions [14, 15, 16, 17, 
18, 19, 20] have been discussed (see also [21, 22]). Following the analogy suggested at the 
beginning of the present script, it is possible to find bag-like solutions in five- dimensional 
geometries with well defined AdSs limit for large values of the bulk radius. It has been shown 
that multi-soliton solutions arise naturally also in a five-dimensional context and fully an- 
alytic solutions have been reported, always in warped geometries, where the field profiles 
describe both kink-antikink systems and trapping bag solutions [23]. 

These discussions have been conducted, initially, in the framework of five- dimensional 
Einstein-Hilbert gravity. Possible extensions of this approach include 

• the modification of the underlying gravitational theory (for instance from Einstein- 
Hilbert gravity to Gauss-Bonnet or Brans-Dicke gravity); 
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• the increase of the dimensionality of the space-time (for instance from kinks in five 
dimensions to Abelian vortices in six dimensions); 

• the increase in the number of defects (for instance from single-defect solutions to multi- 
defect solutions). 

One of the distinctive aspects of the aforementioned attempts is that the obtained field 
profiles are thick, i.e. they have an internal structure and they are obtained as a consistent 
solution of the underlying gravitational theory. In this sense the nature of the physical 
system greatly differs from their (1 + 1) -dimensional counterpart. Particularly interesting 
situations are the ones where the geometry is completely regular 2 and it has well defined 
AdSs limits for large absolute value of the bulk radius. 

The gravity theory employed in the actual constructions plays a relevant role in the 
physical properties of the solutions. As already mentioned, the obvious choice is to embed 
the defects in a higher-dimensional Einstein-Hilbert gravity theory (see, for instance, [20] and 
references therein). However, also different choices are equally plausible. For instance one 
might choose to select a theory where the Einstein-Hilbert terms is non-minimally coupled 
to the scalar degree of freedom [24, 25, 26, 27]. Gravitating kinks have been also studied in 
quadratic gravity theories where the Einstein-Hilbert action is complemented by the Gauss- 
Bonnet combination [28, 29, 30, 31]. It is finally possible to combine these two options and 
discuss the possibility of higher-dimensional solutions where the Einstein-Hilbert action is 
complemented by the Gauss-Bonnet combination and the dilaton couples explicitly to the 
curvature in the string frame metric. This possibility has been developed in [32, 33, 34]. 

Another arena for thick defects involves also the increase of the dimensionality of the 
space-time. In the framework of the Abelian-Higgs model it is possible to find thick string 
solutions in a six-dimensional space-time with well defined AdS 6 limit far from the vortex 
and a Minkowskian limit close to the core of the defect. These solutions may arise, with 
different features, both within Einstein-Hilbert gravity [35, 36, 37, 38, 39] and within Gauss- 
Bonnet gravity [40]. Another example is represented by higher dimensional hedgehogs which 
are present in the case of seven (warped) dimensions (see [41, 42] and references therein). 
Finally, there exist, always in a five-dimensional context, solutions where the solitons have 
non-topological properties [43], i.e. the static fields have a bag-like profile as a function of 
the bulk radius. 

In the present paper the main theme will be to investigate time- dependent soliton so- 
lutions that are inspired by the symmetries of the static solutions. There is the hope, in 
this context, of finding new classes of cosmological models that do not follow in the class 
of separable solutions that have been already exploited in the five-dimensional case. The 
latter solutions are typically obtained in the presence of thin brane configurations that arise 

2 In this paper, when talking about regular geometry, we will always refer to the case when the relevant 
curvature invariants are finite for all values of the coordinates lying in the domain of definition of the 
solutions. 
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as Dirac delta functions in the energy-momentum tensor of the brane. Typically these so- 
lutions are singular at the location of the brane and the defects, being thin, do not have 
internal structure. In different words we can say the the interest of the present exercises 
resides mainly in the fact that time- dependent solitonic solutions can be eventually used to 
describe, for a fixed value of the bulk radius, higher- dimensional cosmological models which 
are, really and truly, stemming from thick gravitating defects or, for short, thick branes. 

To introduce the discussion in more specific terms consider the static solutions of the 
following five- dimensional theory 3 : 



~ + \G**d A tpd B tp-V{tp) 



XI) 



S = J d 5 x^f\G\ 

in the case of the line element 

ds 2 = G AB dx A dx B = a 2 (w)[r]^dx^dx u - dw\ (1.2) 

where w is the bulk coordinate, a(w) is the (static) warp factor and rj^ is the four- 
dimensional Minkowski metric 4 . Consider then, a specific solution of this theory (see 
[14, 16, 18] and references therein), i.e. 

a H = m \ , i ' ^ = V arctan M» (1-3) 

\ , X 2 W I + 1 V K 

arising when the potential V(<f>) takes the following trigonometric form: 



V(<p) = — 



1-5 sin 2 $ 




(1.4) 



From Eqs. (1.3) and (1.4) we can infer two kinds of solutions: 



• boosted solutions that can be obtained by replacing w with / -f(u)(w + ut) or with 
~f(u)(w — ut) (where ^(u) = (1 — w 2 ) -1 / 2 ); 

• time-shifted solutions. 

The first class of solutions describes kinks that travel along the fifth coordinate and they 
will be discussed a bit later in this paper. The second class of solutions can be more directly 
interpreted as a time-dependent cosmological model for a fixed value of the bulk coordinate 
(that can be chosen to be, for instance, at w — 0). The trouble with this second class of 
solutions is that, in general, they will not be diagonal. Indeed, from Eqs. (1.2), (1.3) and 



3 In Eq. (1.1) k = 8ttG5 = 8w/M 3 and G5 is the five-dimensional Newton constant. 

4 Conventions: capital Latin indices run over the whole dimensionality of the space-time while Greek 
indices run over the four-dimensional space-time; Latin (lowercase) indices run over the three-dimensional 
spatial geometry. The signature of the metric is mostly minus, i.e. (+,—,—,—,—) 



3 



(1.4) it is possible to deduce that a time- dependent solution of the theory can be found, 
always for the potential (1.4), for a line element of the type 5 



ds 2 = a 2 (w, t)[cIt 2 — dx 2 — dw 2 — 2dw dr], 



(1.5) 



where 



a(w, t) 



\ 2 (w + t) 2 + 1 



1 




arctan [X(w + r)], 



(1.6) 



The easiest way of getting convinced of this statement is by using the symmetry of the metric 
and then, as a cross-check, to verify that the transformed metric is a solution of the action 
(1.1) with the potential given in Eq. (1.4). Indeed, the Ricci tensor obtained from Eqs. (1.5) 
and (1.6) has covariant components that can be computed to be: 



By contraction of the Einstein equations derived from Eq. (1.1), the relevant system can be 
recast in the following form: 



where is the covariant derivative computed in terms of the five- dimensional metric. Using 
the field profile appearing in Eq. (1.6) and the potential of Eq. (1.4) into Eqs. (1.8) and (1.9) 
the system is satisfied and the solution is proven to exist. Notice that the off-diagonal form 
of the metric reported in Eq. (1.5) is essential to match the off-diagonal components of the 
Ricci tensor (i.e. Rq w and R w o) which are now non-vanishing since the warp factor depends 
explicitly upon time. Indeed, from Eq. (1.5), the Ow and wO components of Eq. (1.8) lead 
also to a generalized momentum constraint that mixes derivatives of the field profile with 

5 In what follows the time coordinate shall be denoted by r. This is motivated by the fact, as it will be 
shown, that the geometries written in the parametrization adopted in the present script reduce to what we 
would call a conformally flat metric for a fixed value of the bulk coordinate. In this case the time coordinate, 
i.e. r will be what is customarily called, in cosmological applications, conformal time coordinate (as opposed 
to the cosmic time coordinate where all the clocks are synchronized). This is the reason why we prefer to 
use r rather than t to denote the time. 



R 




(1.7) 




Rab = « d A ipd B ip - -V , 



(1.8) 



(1.9) 
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respect to w and with respect to r. The momentum constraint is solved iff the metric is in 
the form of Eq. (1.5), i.e. off-diagonal. 

This example is, at the same time, tangential and inspiring. It is tangential since, as 
already remarked, the obtained geometry is necessarily off-diagonal and, therefore, difficult 
to interpret as a cosmological metric at fixed value of the bulk coordinate. It is, however, 
also inspiring since it suggests that it might not be hopeless to find time-dependent soliton 
solutions leading to a non-separable form of the warp factor. By separable form of the warp 
factor we mean solutions where a(w, r) can be written as the product of two functions one 
depending solely upon w and the other depending solely upon r. Separable solutions are 
normally discussed in the framework of thin brane models but, as we will argue, they are 
rather unnatural in the case of thick defects. 

The plan of the forthcoming discussion will therefore be the following. In section 2 we 
will look for diagonal (time-dependent and non separable) solutions in the five- dimensional 
case and with a single scalar profile. In section 3 this discussion will be generalized to the 
case of mult i- defects. Section 4 contains some complementary material on traveling solitons. 
Finally section 5 contains our concluding remarks. 



2 Time-dependent gravitating defects 



To improve on the off-diagonal ansatz shown to satisfy the whole system of equations sup- 
ported by time-dependent (gravitating) profiles discussed in the introduction, it is appropri- 
ate to consider the following ansatz: 



ds 2 = G AB dx A dx B = a 2 (w,r)[dT 2 - dx 2 } - b 2 (w,r)dw 2 . 



(2.1) 



Using Eq. (2.1) for the calculation of the Ricci tensor and of the Christoffel connections the 
various components of Eq. (1.8) easily become, after some algebra: 



— [W + AH 2 - HT\ -[F + 3H + F(F - H)] = k ^ 
tr L 

a 2 ■ 2 

-tt[W' + 4H 2 - HT\ + [H + 2H 2 + HF] = -nVa 2 , 
fr 3 



-W 



: [F + F 2 + 2HF] - AH' + m{T -H) = k 



V ' 2 + ^Vb 2 



H-4H' + 3FH = ivjxp'. 
The notations employed in Eqs. (2.2), (2.3), (2.4) and (2.5) are the following: 

• the prime denotes a derivation with respect to the bulk radius w; 

• the over-dot denotes a derivation with respect to the conformal time coordinate r; 



(2.2) 

(2.3) 

(2.4) 
(2.5) 
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• 7i and T are related with the derivatives of a and b with respect to w, i.e. more 
specifically, 

ni f x <91na 9 In 6 

^,r) = — , %r)^; (2.6) 

• /J and F are related with the derivatives of a and b with respect to r, i.e. more 
specifically, 

d In a ri . . (9 In b 

H(w, r) = — , F(w, r) = ^r. (2.7) 

Equations (2.2), (2.3) and (2.4) follow, respectively, from the (00), (i = j) and (ww) com- 
ponents of Eq. (1.8) while Eq. (2.5) follows from the (Ow) component of Eq. (1.8). Finally, 
Eq. (1.9) leads, within the notations expressed by Eqs. (2.6) and (2.7), to the following 
explicit relation 

a 2 a 2 (W 

(p ~ ^" + (2H + F)<p - ^(4W - ?y + j^a 2 = 0. (2.8) 

If the dependence upon r and w would be separable (i.e. a(w, r) = a(r)a(w) and b(w, r) = 
P(r)b(w)), then, necessarily, 7i and JF will only be function of the bulk radius (i.e. 7i(w) and 
J-(w)), while H and F will only be functions of the conformal time coordinate (i.e. H{t) 
and F(t)). The constraint of Eq. (2.5) can then be trivially satisfied by requiring either 
F = (i.e. f3 constant in r) or TL = (i.e. a constant in w). This situation is, however, not 
generic. 

It is appropriate to mention that by taking linear combinations of some of the above 
equations, the resulting expressions can be put in a more friendly form. In particular, by 
summing up Eqs. (2.2) and (2.3) we do get: 

(p 2 = -[2H 2 + 2HF - F 2 -F-2H}. (2.9) 

K 

By combining Eqs. (2.3) and (2.4) we do get instead 

p' 2 = -(HT - W) + -^r[(F -H) + (F + 2H)(F - H)}. (2.10) 

K Kb 

Consider then the following form of the warp factors: 

a(w,r)= 1 =, b(w,r)= t =, (2.11) 

tXHw + rY + l J\ 2 (w + t) 2 + 1 



where e is both constant in w and in r. It can be verified that Eq. (2.11) satisfies the whole 
system of equations provided the scalar field and the potential are: 

(p(w,T) = \-aictm[X(w + T)], (2.12) 

V K 

3A 2 e 2 - 1 . „ T , x , [k 



V(<P) = -t ^(5sin 2 <l»-l), $ = J^. (2.13) 



2k e 2 
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Figure 1: The warp factors for different values of the integer \i. 



If e — > 1 the solution does not reduce to a static limit but rather it reduces to a free solution 
with vanishing potential (i.e. V(ip) = 0) where the explicit expression of ip is still given 
by Eq. (2.12). Again, the solution e = 1 does not have a static analog since the time 
dependence is essential. In other words $(«;) = arctanAw is not a solution of the system if 
we set the time dependence to zero and V(<&) = from the very beginning. The solution of 
Eqs. (2.11), (2.12) and (2.13) can be generalized to a family of solutions characterized by 
an integer parameter \x. The solution can be written as 

a(w,r) = — , b(w, r) = ea(w, r), (2-14) 

[A 2 ^(w + r) 2 ^ + 

where, as before, e is a constant. The form of the warp factor given in Eq. (2.14) solves the 
whole system provided 



(p{w,r) 
V(<p) = 



'3(2/1 



K 



■ arctan [A M (w + r) 



3A 2 /e 2 



2k 



where \i > 1 and 



AH 



sin 2 $, 



(2.15) 
(2.16) 

(2.17) 



/ 3(2^ - 1) 

In the case \i = 1 the solution expressed by Eqs. (2.14), (2.15) and (2.16) goes exactly to the 
solution of Eqs. (2.11), (2.12) and (2.13). Again, in the limit e = 1 Eqs. (2.14) and (2.15) 
still solve the system with vanishing potential. 

Different values of fi imply different properties for the topology of the field solution. 
This point is illustrated in Figs. 1, 2 and 3 where the warp factors and the field profiles 
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Figure 2: The field profile for odd values of ji. Natural gravitational units 2k = 1 have been 
adopted. 

are reported for different values of the integer /i that appears in Eqs. (2.14), (2.15) and 
(2.16). According to Fig. 1, when fi increases from 1 to 4 the warp factor becomes broader 
in the (w,t) plane (recall that b(w,r) = ea(w,r)). For fixed value of the conformal time 
coordinate, say r = 0, a(w,0) — > |Aw| _1 for \Xw\ — > oo and b(w,0) — > e|Aw| _1 . In this case 
the space-time geometry can be viewed, for large absolute value of the bulk radius, as a 
slight deformation of AdSs. The deformation parameter is exactly e, which is dimensionless. 
As soon as e — ► 1 the exact AdS^ limit is recovered for \w\ — > oo. The geometry is always 
regular and all the curvature invariants do not have any pole for any finite values of either w 
or r. In the limit e — > 1 the Weyl invariant vanishes. A change in the value of fi does affect 
crucially the topological nature of the field profile. This aspect can be simply understood 
from Eq. (2.15). If fi is even the field profile interpolates between two different minima of 
the underlying potential and it has the characteristic kink-like shape. In Fig. 2 this aspect is 
illustrated for two different odd values of fi, i.e. fi = 1 and \x = 3. Notice a slight difference 
between the case /i = 1 and the case fi = 3. In the case /i = 3 (see Fig. 2, plot at the right) 
the profile flattens around the origin of the coordinate system. In fact, in the case \x = 1 the 
first derivative of r) (both with respect to w and r) is not vanishing at the origin, i.e. 
for w = t = 0. On the contrary, at the same point, in the case = 3 the first and second 
derivatives of <&(«;, r) are both equal to and the only non- vanishing derivative at the origin 
will be the third one. Besides this point we can safely say that, when fi is odd, the field 
profile exhibits a topological nature. Conversely, when fi is even, the field profile does not 
have a kink-like structure but rather a bag-like structure signaling that the corresponding 
field profile is non-topological. This aspect is usefully illustrated in Fig. 3 where (p(w,r) is 
reported in the cases \x = 2 and // = 4. As it appears from Eq. (2.15), when /i is even, the 
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Figure 3: The field profiles for even values of ji are reported in the natural gravitational 
units introduced in Fig. 2. 

field profile tends to the same constant values when w and r are large in absolute values. 
As far as the cosmological properties of the solutions are concerned, it is useful to look at 
Fig. 4 where H and H are reported as a function of the conformal time coordinate r. It is 
clear from Fig. 4 that, for fixed value of the bulk radius, say w — 0, H(t) changes sign once 
and H changes sign twice. This behaviour is typical of a well known class of cosmological 
scenarios called bounces. Bouncing solutions in the Einstein frame (which is the one adopted 
in the present analysis) are solutions where the Hubble rate flips its sign [44, 45] (see also 
[46] and references therein). The bouncing behaviour obtained in the present solutions and 
illustrated in Fig. 4 is, at the same time, unrealistic and encouraging. It is unrealistic since 
a physical bounce should have features that are opposite to the ones exhibited by Fig. 4. In 
other words a realistic bounce should have H < for r < and H > for r > 0. This type 
of evolution would imply that the contraction took place in the far past (for r < 0) while 
today, as we know, the Universe expands rather than contracting. In Fig. 4 the behaviour 
of H (and H) is opposite: the Universe expands in the past and contracts in the future, 
which is unrealistic. At the same time, however, the obtained results are encouraging since 
we showed, analytically, that a bouncing behaviour arises naturally. 

It is relevant to mention that the only separable solutions that are allowed by (single) 
gravitating solitons are just the trivial ones where H = and F = 0. To demonstrate this 
point consider the following ansatz for the warp factors: 

aM = 7m- bM= 7m- (2 - 18) 

supplemented by the condition <p = <p(w). Using Eq. (2.18) into Eqs. (2.6) and (2.7) it 
follows that Ti = and H' = 0. But then the momentum constraint of Eq. (2.5) implies 
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Figure 4: The conformal time evolution of H and H is illustrated at a fixed value of the 
bulk radius and for different values of fi. 



that, since = 0, also F — 0. The requirement that the solution is separable demands, 
moreover, that all the pieces containing H, F and their derivatives vanish in Eqs. (2.2), 
(2.3), (2.4) and (2.8). This condition translates, after some algebra, in the following pair of 
equations 

H + H 2 = 0, HF + H 2 = 0. (2.19) 

But since F = (because of the momentum constraint) it must also happen, according to 
Eq. (2.19) that H = 0. Thus, a(r) and f3{r) can only be constants. 



3 Time-dependent gravitating multi-defects 



The considerations developed in the previous section can be generalized to the situation 
where two scalar degrees of freedom are simultaneously present. The action to be discussed 
in the present section will therefore be: 



S = J d 5 x^\G\ 



~ + l -G AB d A ^d BV + l -G AB d AX d B x - W(<p, x ) 

2k l 2 



(3.1) 



Using the same conventions of the previous section the evolution equations only involving 
<f, x an d their first derivatives can be written, in the geometry defined by Eq. (1.5), as 



y? 2 + x 2 = -\2H(H + F)-2H-F 2 - F], 

K 

W' + XX' = -[W-4:H' + 3FH}, 

K 

b 2 



/2 , i2 



¥ +x 



K 



\HT -rC) 



[(F — H) + F + HF — 2H 2 }. 



(3.2) 
(3.3) 
(3.4) 
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The remaining three equations involve directly the scalar potential W(tp, x) and they are: 



W(<p, X ) 

,2 



2na 2 



[H + 2H 2 + HF) 



2Kb 2 



(H' + m 2 - HF), 



(P ~ + [(2H + F)<p - a - 2 {m - FW] + a 2 ^ 



= 0, 



a 2 ,, x a 2 , , ., cW 

X - ^x" + [(2/T + F) x - ¥ W- F) X '] + a 2 — = 0. 



(3.5) 
(3.6) 
(3.7) 



Consider, therefore, the following ansatz for the field profiles and for the warp factors: 



tp(w,T) = vyj\ + g(w,r), x(w,t) = vyjl - g(w, r), 

b(w, r) = ea(w, r). 



a(w,T) 



1 - g 2 {w,T) 



(3.8) 
(3.9) 



where g(w,r) is a monotonic function which is continuous with its first and second derivatives 
with respect to w and r. Furthermore, we do assume that \g(w,T)\ < 1 so that the square 
roots that appear ubiquitously in the ansatz are always well defined in the domain where w 
and r are allowed to vary. By inserting Eqs. (3.8) and (3.9) into Eqs. (3.2), (3.3) and (3.4) 
the following conditions on g(w, r) can be explicitly obtained: 



9 = 



3ggg' 



9 



„9 2 
9 = 9^. 



(3.10) 



To obtain Eq. (3.10) the arbitrary constant v has been expressed in terms of k as kV 2 = 6. 
The solution of Eq. (3.10) becomes then: 



X(w + r) 



X 2 (w + r) 2 + l 



(3.11) 



Equations (3.5), (3.6) and (3.7) can then be used to determine the scalar potential that 
becomes, within our conventions, 



W (e,v 



^V + x 2 )-(^ 4 + x 4 ) 



+ W 1 (e,v)( l f 2 + X 2 -2V 2 ) 



y(y? 4 + x 4 )-^V + x 2 ) 



where 



W (e,v) 



15\ 2 /l-t 



v 4 \ 



W 1 (e,v) 



2vA 



(3.12) 



(3.13) 



It can be directly checked that Eqs. (3.12) and (3.13) solve directly Eqs. (3.5), (3.6) and 
(3.7). The physical situation described by the reported solution is the one where we have a 
kink and an anti-kink both time dependent. This aspect is illustrated by Fig. 5 where the 
functions (p(w, r) and x( w i T ) are separately reported in the plot at the left and in the plot 
at the right. 
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Figure 5: The kink-antikink system of Eq. (3.8) in the case v — 1 and always in natural 
gravitational units. 

As it was shown in the previous section, the one-field system only allows separable solu- 
tions which are trivial. If two fields are present the situation is different. Consider, for in- 
stance, the situation where the warp factors of the line element of Eq. (1.5) are parametrized 
as 

= 73?6tT' (3 - 14) 

namely b(w) is just static while the time dependence of a(w, r) is separable. In this situation 
Eq. (3.3) is trivially satisfied while Eqs. (3.2) and (3.4) lead, respectively, to the following 
pair of conditions: 



y/ 2 = \<H 2 - W) - -^[H + 2H% (3.15) 
X 2 =-(H 2 -H). (3.16) 

K 

The potential can then be written as 

Wfa X ) = ^ 2 (H + - ^(W + m% (3.17) 

There are now different possibilities and, for purposes of illustration, let us focus the attention 
on the simplest one, namely the case where H + 2H 2 = 0. In this situation the potential 
will be effectively independent on x an d the evolution of at(r) will be simply driven by the 
kinetic energy of x- in i&ct, Eqs. (3.15) and (3.16) reduce to 

V K A Z W Z + 1 

X=f-H{r) (3.19) 

V K 
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and H(t) = (2r) _1 which is, as expected, the expansion rate dictated typical of a FRW 
Universe filled by a free (minimally coupled) scalar degree of freedom. Since H = —2H 2 , 
Eq. (3.17) allows to determine also the potential of (p. The scalar profile <p(w) and its related 
potential coincide, therefore, with the ones reported, respectively, in Eqs. (1.3) and (1.4). 



4 Traveling solitons 

In this section we are going to discuss the case of traveling soliton solutions. More specifically 
it will be shown that, under certain conditions, every static solution cane be promoted, by 
Lorentz invariance, to traveling soliton solution. To discuss this point it is useful to write 
some of the equations of the previous section in terms of the new boosted coordinates: 

w + Ut w -Ut 

U = , V = ; (4.1) 

Vi- u 2 Vi -u 2 v ; 

Consider therefore the line element of Eq. (1.5) with a(w,r) = b(w,r). After few linear 
combinations it is easy to show, with simple algebra, that Eqs. (3.2), (3.3) and (3.4) translate 
in the following triplet of equations: 

<fil + xl = l(hl-h uu ), (4.2) 
ti + xl = -(hl-h vv ), (4.3) 

K 

3 

<Pu<Pv + XuXv = -{KK - h uv ), (4.4) 

K 

where the subscript denotes a partial derivation with respect to the corresponding variables; 
moreover h — In a. Equation (3.5) implies, in terms of the variables o Eq. (4.1), 

3 

W(<p,x) = -—[(h uu + h vv + 3h 2 u + 3h 2 v ) + 2cosh2a(h uv + 3h u h v )}, (4.5) 
2k 

where cosh 2a = (1 + U 2 )/(l — U 2 ). It is finally easy to show that Eqs. (3.6) and (3.7) lead, 
respectively, to the following pair of equations: 

2 dW 

Vuv - ¥uu - ¥w - 3(h u (p u + h v ip v ) - 3 cosh 2a(h u ip v + h v (p u ) + a — — = 0, (4.6) 

ocp 

Xuv - Xuu - Xw - KKXu + KXv) - 3 cosh 2a(h u Xv + KXu) + a<2 ~^ = °> ( 4 - 7 ) 

Suppose now that both, (p and \ are functions of u, i.e. = ip(u) and \ = Xu- Then, Eqs. 
(4.3) and (4.4) will be satisfied provided h = h(u), i.e. a = a(u). Under these conditions the 
relevant set of equations to be solved becomes: 

vl + xl = - K {h 2 u -h uu ), (4.8) 
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W(ip, X ) = ~(h uu + 3hl) (4.9) 
Vuu + 3h u ip u - = 0, (4.10) 

Xuu + ZKxu - a 2 ^- = 0. (4.11) 

It can be easily seen that Eqs. (4.8), (4.9), (4.10) and (4.11) have the same form of the 
static solutions (in the two field case) when a(w,r) = b(w,r). The derivative with respect 
to w translates into a derivative with respect to u. So, every static solution can be boosted 
and the resulting metric and field profiles will still be a solution. The same conclusion can 
be reached if the boost is performed in the opposite direction. In this second case the field 
variables will only be functions of v (and not of u) and the relevant equations can be obtained 
from Eqs. (4.8)-(4.11) by replacing u — > v in the derivatives of 0, x an d h. 

As an example it is useful to consider the following solution of the above field equations 
in the case where, both field profiles are only functions of, say, u: 

= W,{[l + 9(«)] S/2 + [l-?(«)] 3/2 }, 
X (u) = Xo {[l + ?(«)] s/2 - [1 - ?(«)] 3/2 }, (4.12) 

where q(u) is a continuous and differentiable function of u; ipo and Xo are both constants. 
As the function g(w, r) introduced in the previous section, it will be required that \q(u)\ < 1 
in the whole domain of definition. To satisfy the constraint of Eq. (4.4) it will be necessary 
that also a = a{u) and the static solution implies that the appropriate ansatz is nothing but 

a(u) = -. (4.13) 

(\ 2 ^U 2 ^ + 1) 2 M 



The function q(u) is therefore determined to be 

q(u) = ^ arctan [W], Kcp 2 = (4.14) 

with xo = V^o- Equations (4.9), (4.10) and (4.11) allow to determine the potential: 

W( V , X ) = ^[sin 2 a] V [(2/i - 1) - (2/* + 3) sin 2 v] + [\<p + x| 2/3 + \<f>- x\ 2/3 ~ 2]A, (4.15) 
where, as usual, 

varphi arphi = 2(p {p and x — ^VoX- The functions a(<f, x) and A(<p, x) appearing in Eq. 
(4.15) are defined, respectively, as: 



a — 2 



A(ip, x) = P^-cos(7{(j[sinV] ^ [(fi- 1) - (2/i + 3) sin 2 a 



7T 



+ ^coso-[sin 2 cr] ? }, (4.16) 
^X) = - A \W + X\ 2,3 -W-X\ 2I % (4.17) 
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As expected this traveling solution has a static analog that has been discussed in [43]. 

To conclude the present section let us consider the possibility of having solutions that 
mix u and w. There are, of course, trivial solutions to this problem, namely the ones where 
the dependence of the various fields and of the warp factors depends either upon {u + v) 
or upon (u — v). In these two cases the resulting profiles are either completely static or 
completely time-dependent. 

The kind of mixed solutions which would be interesting to obtain would be the fivc- 
diemensional analog of the (1 + 1) dimensional solutions of the Sine-Gordon system that 
have been introduced in [1, 2] and extensively discussed in [3]. Indeed, it is well known 
that there are no static two-kink solutions in the Sine-Gordon system when only one scalar 
degree of freedom is present. However, two-kink solutions in which either kinks have arbitrary 
velocity (but at least one non-zero) can be constructed [1, 2, 3]. A typical example of this 
dynamics is the field profile given, in (1 + 1) dimensions, by 4 arctan [U sinh (72;) / cosh (■yU t)} 
where 7 = (1 — U 2 )^ 1 ^ 2 and the two relevant coordinates have been denoted by x (space) and 
t (time). In this solution the kink approach the x axis with velocity U, scatter elastically 
at t — and emerge from the interaction with a computable phase shift. Notice also that, 
for t — > +00 and x — > +00 the field profile becomes solely function of u — 7(2; + Ut). 
Conversely, for x — > +00 and t — > —00, the field profile becomes, asymptotically, only 
function of v — 7(2; — Ut). 

In the five-dimensional warped geometries discussed in the present paper these solutions 
are rather hard to find and we cannot report any positive result. The main problem is 
represented by the structure of Eqs. (4.2), (4.3) and (4.4). Indeed it seems rather difficult to 
satisfy simultaneously the momentum constraint of Eq. (4.4) together with the conditions 
implied by Eqs. (4.2) and (4.3). In spite of this negative result it seems to be worthwhile to 
investigate further the possible solutions describing the elastic scattering of solitons in five 
dimensions. These solutions might have an intriguing cosmological interpretation. 

5 Concluding remarks 

In the present paper time-dependent solitonic solutions have been discussed. It has been 
shown that these solutions exist both for the case of topological solitons and for the case of 
non-topological solitons. The discussion has been also generalized to include the possibility 
of multi-defects, i.e. the situation where more than one defect is present in the system. In 
many respects the results of the present analysis are encouraging but still preliminary. They 
are encouraging since there exist solitonic profiles leading to time-dependent geometries that 
are diagonal and that have a well defined AdS 5 limit for large values of the bulk radius and 
for a fixed value of the conformal time coordinate. The reported results suggests several 
lines of possible developments. In the first place more realistic bouncing behaviours may 
be studied always in the case of five-dimensional warped geometries. Moreover, it is not 
excluded that different classes of solutions (like the ones related to elastic scattering of two 
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profiles) may emerge. 

The main aim of the present script has been to provide examples of thick defects that 
lead to solutions whose time-dependent evolution may be potentially relevant for cosmology. 
In this respect, an obvious extension of the reported considerations will be to study higher 
dimensional defects, like, for instance, Abelian vortices in six-dimensions. These models are 
known to provide, at the static level, AdS 6 geometries for large values of the bulk radius. 

It would be finally interesting to compute the production of relic gravitons. Indeed, 
while at the static level the zero modes of the gravitons will certainly be localized owing 
to the AdSs nature of the geometry, if the solutions exhibit time-dependence the natural 
question to ask will be how many gravitons are produced by the time-dependent action of 
the gravitational field. This analysis requires necessarily a more realistic construction: it 
would not be enough, in the present approach, to patch different solutions since a more 
realistic implementation of the bouncing behaviour is required. All these issues are beyond 
the aims of the present investigation. 
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